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ABSTRACT 

We present a comprehensive study of the star cluster population and the hierarchical structure in the clustering 
of blue stars with ages ^ 500 Myr in the Local Group dwarf irregular galaxy NGC 6822. Our observational 
material comprises the most complete optical stellar catalog of the galaxy from imaging with the Suprime-Cam 
at the 8.2-m SUBARU Telescope. We identify 47 distinct star clusters with the application of the nearest- 
neighbor density method to this catalog for a detection threshold of 3cr above the average stellar density. The 
size distribution of the detected clusters can be very well approximated by a Gaussian with a peak at ~ 68 pc. 
The total stellar masses of the clusters are estimated by extrapolating the cumulative observed stellar mass 
function of all clusters to be in the range 10 3 - 10 4 M©. Their number distribution is fitted very well by a 
power-law with index a ~ 1.5 ± 0.7, which is consistent with the cluster mass functions of other Local Group 
galaxies and the cluster initial mass function. 

The application of the nearest-neighbor density method for various density thresholds, other than 3er, enabled 
the identification of stellar concentrations in various length-scales, in addition to the detected star clusters of 
the galaxy. The stellar density maps constructed with this technique provide a direct proof of hierarchically 
structured stellar concentrations in NGC 6822, in the sense that smaller dense stellar concentrations are located 
inside larger and looser ones. We illustrate this hierarchy by the so-called dendrogram, or structure tree of the 
detected stellar structures, which demonstrates that most of the detected structures split up into several sub- 
structures over at least three levels. We quantify the hierarchy of these structures with the use of the minimum 
spanning tree method. We find that structures detected at 1, 2, and 3a density thresholds are hierarchically 
constructed with a fractal dimension of D « 1.8. Some of the larger stellar concentrations, particularly in the 
northern part of the central star-forming portion of the galaxy, coincide with IR-bright complexes previously 
identified with Spitzer and associated with high column density neutral gas, indicating structures that currently 
form stars. The morphological hierarchy in stellar clustering, which we observe in NGC 6822 resembles that of 
the turbulent interstellar matter, suggesting that turbulence on pc- and kpc-scales has been probably the major 
agent that regulated clustered star formation in NGC 6822. 

Subject headings: galaxies: dwarf — galaxies: individual (NGC 6822) — galaxies: irregular — galaxies: star 
clusters — open clusters and associations: general 



1. INTRODUCTION 

Star formation in galaxies is observed on a wide range 
of scales, from stellar complexes and aggregates over OB 
associations to compact embedded clusters, which again 
often show substructure. These systems are not distinct, 
independent entities but rather appear to form a contin- 
uous hierarchy of structures over all these sca les (e.g. 
Efremov & Elmegreen 1998; Elmegreen et al. 2000). The in- 
terstellar matter (ISM) also shows a hierarchical structure 
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from the largest giant molecular clouds down to individual 
clumps and cores. The complex hierarchical structure of the 
ISM is believed to be shaped b y supersonic turbulence (e.g., 
Ballesteros -Paredes et al.1 120071). The sca ling relations ob- 
served in molecular clouds (lLarsonlfT98rl) can be explained 
by the effect of turbulence, where energy is injected at very 
large scales and cascades down to the smallest scales, creat- 
ing eddies and leading to a clumpy, filamentary structure on 
all scales. These structures appear self-similar and are there- 
fore often described as fractal. Turbulence is believed to play 
a major role in star formation by creating density enhance- 
ments that become gravitationally unstable and collapse to 
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form stars. The spatial distribution of young stars and star 
clusters probably reflects this process. 

There are only few cases, where the formation of stellar 
structures in star-forming galaxies with resolvable stars can 
be studied in the complete extent of the galaxy. NGC 6822 
is one of these rare galaxies, because it is far enough from 
us for its whole extent to be fairly covered in small fields- 
of-view, but also close enough for its bright stellar content 
to be sufficiently resolved. NGC 682 2 (IC 4895) is a Loca l 
Group dwarf irregular of type Ir IV-V (Ivan den Berghl feOOO). 
It is the closest dwarf irregular apart from the Magellanic 
Cl ouds (MCs), located at a distance of 490 ± 40 kpc from 
us dMateoll 19981) . and belongs to an extended cloud of irreg- 
ulars named the "Local Group Cloud". Due to its close dis- 
tance NGC 6822 appears quite extended on the sky; its opti- 
cal angular diameter is over a quar ter of a degree, while it s 
Hi disk measures close to a degree (Ide Blok & Wa lter 2000). 
NGC 6822 is relatively near the Galactic plane (b = -18.4°) 
and therefore suffers from si gnificant foregrou nd extinction 
from the Milky Way (MW). lMassev et all (120071) using multi- 
band imaging of NGC 6822 within the Survey of Local Group 
Galaxies Currently Forming Stars, derived a total interstellar 
reddening of E(B-V) = 0.25 ±0.02, with E(B-V) = 0.22 be- 
ing Galactic. Other recent optical studies of the stellar conten t 
of the galaxy dde Blok & Walter! 12003; Batti nelli et aT1l2003l) 
show that NGC 6822 has an extended stellar distribution, with 
the young blue stars foll owing the distribution of the Hi disk 
( Komiyam a et al. ll2.Q03l> . while the old and intermediate-age 
stellar popu lation is found significa ntly more extended than 
the Hi disk (de Blo k & Walterj|2006l) . 

The stellar and gas components of NGC 6822 could not 
explain the shape of the high-resolution rotation curve ob- 
tained with the ATC A, except of the very inner regions 
(IWeldrake et al.l 120031) . and therefore it is considered to be 
very dark-matter-dominated. NGC 6822 is a metal-poor, 
relatively gas-rich galaxy with ISM metal abundance of 
about 0.2 Z(T), i.e. the t otal fraction of metals is Z ~ 
0.004 (Skillma iTetaI1ll989j) . and total Hi mass of - 1.3 x 
10 8 M g dde Blok & Walter! 120001) . Its star formation rate 
(SFR), based on Ha a nd far-IR fluxes, is found to be 
around ~ 0.06 M yr"' dMateol 119981; Hsrael et alj [7996b or 
^0.01 M Q yr _1 dHunter & Elm egreen 2004). Evidence for in- 
cr eased star form ation be tween 75 and 100 M yr ago is found 
bv iHodgdd 19801) . while iGallart et alj dl996al) found that the 
star formation in NGC 6822 increased by a factor of 2 to 6 be- 
tween 100 and 200 Myr ago. Moreover, HST imaging showed 
that the rece nt SFR is spat ially variable in the central parts of 
NGC 6822 dWyderl F2001). All these findings are consistent 
with the mostly constant but stochastic recent star formation 
histori es often derive d for other dwarf irregular galaxies (see, 
e.g. JWeisz etai1l2008l) . 

NGC 6822, like most dwarf irregulars, provides an ideal 
laboratory for the study of galaxy evolution, and in particular 
at early stages, since its low metallicity and rich gas content 
suggest that NGC 6822 is still in an early stage of its con- 
version from gas into stars. Moreover, its relatively simple 
structure, without dominant spiral arms or bulge, makes the 
study of the various physical processes related to star forma- 
tion relatively straightforward. The scope of the present study 
is the comprehensive understanding of the structural behavior 
of star formation during the period of the last ~ 500 Myr over 
the whole extent of NGC 6822. Specifically, we investigate 
the cluster population of the galaxy in terms of the size, mass 
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FIG. 1.— Three views of NGC 6822. Top: Integrated ATCA Hi column 
density or zeroth-moment map of NGC 6822. The contour indicates the edge 
of the Hi disk of the galaxy at 5 X 10 20 cirT 2 . The grayscale levels run up to 
4.3 X 10 2 ' cirT 2 (black), which is the maximum column density occurring in 
this map. The beam of 42"4 X 12"0 is indicated in the bottom right. This 
col umn density map is ad apted from dBW06, for additional representations 
see Weldrake et al. 1 2003). Middle: Spatial distribution of the blue stars with 
B—R < 0.5 on the sky. This distribution highlights the concentration of the 
blue stars with ages $ 500 Myr in the inner parts of the galaxy, and the 
general extent of these stars in a disk-like structure, similar to the Hi disk 
of the galaxy. Bottom: Number density map of the blue stars shown in the 
middle panel, measured in 64" X 64" boxes. The iso-density contour levels 
run from the mean background density in steps of la, a being the standard 
deviation of the background density. Isopleths of density > 5<r are plotted 
with thicker lines. 

and structural characteristics of the detected stellar systems. 
We also study the clustering behavior of stars and its relation 
to hierarchically structured stellar concentrations within the 
whole extend of NGC 6822. 
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A comprehensive study of the star complexes, i.e., large- 
scale star-forming regions with size s of the order of few hun - 
dreds pc, was recently presented by Karampelas et al. (2009). 
These authors provide a complete catalog of star complexes 
in NGC 6822, their positions and sizes, based on the data 
collected within the Survey of L ocal Group Galaxies Cur - 
rently Forming Stars (SLGGCFS; iMassev et ai1l2006l 120071) . 
This study is limited to the brightest stars of the galaxy with 
B ^ 22.5 due to the detection limit of the SLGGCFS and 
only to its central part, which is covered by the survey. In 
our investigation here we extend the study of the spatial dis- 
tribution of young stars and hierarchical star formation in 
NGC 6822 to the smallest detectable structures, and to the 
complete extent of the galaxy, including the potentially tidal 
systems, with the use of original photometric material that 
comprises a stellar sam ple 6 times larger than that used by 
Karampela s et al.l ([2009). Our catalog extends to fainter mag- 
nitudes by B ~ 2.5 mag on the main sequence, and covers the 
whole area of the Hi disk of the galaxy (Fig-Stop). As a con- 
sequence, in the present study we investigate star formation in 
a wider range of lengthscales, covering structures with sizes 
from few tens pc (clusters and small associations) to few hun- 
dreds pc (stellar aggregates and complexes). Apart from the 
study of hierarchy, our detection method allows us to also ad- 
dress important issues of star formation, such as the structural 
behavior of the clusters in NGC 6822 and the cluster mass 
function of the galaxy. 

The observational material, which inclu des the most com- 
plete stellar sample observed in NGC 6822 dde Blok & Waited 
2006, from here on dBW06) is described in §|2j In §[3] we de- 
fine the stellar content of the galaxy for the purpose of our 
study and select the type of stars our analysis focuses on. We 
identify the clusters of the galaxy with the application of the 
nearest neighbor (from here on NN) method in § |4j and we 
discuss their sizes and masses in § [5] In § [6] we construct 
the cluster mass function of NGC 6822, and compare our re- 
sults with those of other galaxies. The structural parameters 
of the detected clusters are measured and cross-coiTelated in 
§ [7] The hierarchical behavior in the spatial distribution of 
the blue stellar content of NGC 6822 is investigated in § [8] 
with the detection of stellar structures in various length-scales 
larger than typical clusters, and the construction of their den- 
drogram, to determine the degree of hierarchical clustering. 
We discuss our results concerning hierarchy in § [9] Finally, 
we summarize the findings of this study in § [T0l 

2. OBSERVATIONAL MATERIAL 

2.1. Photometric Stellar Catalog 

In this investigation we use the photometric data derived 
from images of NGC 6822 taken with Suprime-Cam on the 
8.2-m Subaru Telescope at Mauna Kea, Hawaii. Suprime- 
Cam consists of 5 x 2 CCDs of 2048 x 4096 pixels each, 
prov iding a total field of v iew of 34' x 27' with a 0."2 pixel 
size dMivazaki et alj|2002l) . We make use of the results based 
on two sets of archival Suprime-Cam images observed in B, 
R and /. The first set consists of two deep pointings, cover- 
ing the whole of the Hi disk of NGC 6822 (see Fig. Q] top). 
It was observed o n 2001 October 15 and 19, and it is de- 
scribed in detail by Komivam a et al. I ((2003). The second set 
of images was taken in 2000 June, during the early days of 
Suprime-Cam, and it consists of a single shallow pointing to- 
ward the optical center of NGC 6822. dBW06 stress the im- 
portance of the combination of the deep photometry with the 



shallow for a complete understanding of the stellar population 
of NGC 6822, due to the high light gathering power of SUB- 
ARU, which saturates stars down to quite faint magnitudes. 
The present study is based on the comprehensive photometric 
analysis of both data sets performed earlier in dBW06 for the 
investigation of the stellar content of NGC 6822. As a con- 
sequence, much of the information provided in this section 
comes from the results of these authors. We repeat it here for 
reasons of completion. 

2.2. Color-Magnitude Diagram 

In the present study we utilize the final merged photomet- 
ric catalog of stellar sources derived in dBW06 by combin- 
ing the catalogs of stars detected in both deep and shallow 
exposures. This catalog contains 250 237 stellar sources 
identified in all three B, R, I wavelengths. In our analysis 
we limit ourselves to only those "high-quality" objects with 
a photometric uncertainty a < 0.1 mag in all three bands, 
constraining the stellar catalog to 208 894 members. The 
B, B-R color-magnitude diagram (CMD) of these stars is 
shown in Fig. [2] This CMD has a larger d ynamic range than 
those derived from previous investigations (iKomivama et al. 1 
l2003tlBattinelli et alJl200llde Blok & Walteril2003l) . and that 
of SLGGCFS (IMassev et al.ll2006L 120071) . probing to fainter 
magnitudes. By adding in the shallow data our stellar sample 
also covers the O- and B-type star regime of NGC 6822, and 
consequently it is the ideal stellar sample for the investigation 
of the most recent clustered star formation in this galaxy. 

3. YOUNG STELLAR POPULATIONS IN NGC 6822 

The CMD of Fig. [2] contains a mixture of different stellar 
populations. A complete discussion on all observed stellar 
types in NGC 6822 is given by dBW06. In short this CMD 
shows three important distinct components: (i) The vertical 
"blue plume" centered onB-ft ~ 0.3, which consists mostly 
of young stars in NGC 6822. (ii) The red population of stars 
covering the area of the CMD with 1 ^ B-R ^ 1.75 and 
m B <; 23 (named the "red-tangle" by iGallart et alJll996bh . 
which contains mainly the old and intermediate-age stellar 
content of NGC 6822 with ages between 1 and 10 Gyr. (iii) A 
vertical band around B-R ~ 1.5, which consists of contami - 
nant Galactic foreground stars (see also Ma ssey et alJl2007l) . 
Also clearly identified as being due to foreground stars is the 
region with 2 < B-R < 3 and m B > 24.5 (dBW06). 

For our study we focus on the blue populations, and there- 
fore we select our stellar sample from the "blue plume". This 
area in the CMD with B-R < 0.75 consists mostly of young 
main-sequence (MS) and more evolved helium-burning blue 
loop (BL) stars down to an age of roughly 0.5 Gyr. In the 
absence of variable foreground extinction as found toward 
NGC 6822, MS stars would occupy the blue side of the plume 
and BL stars the red. Therefore, since we are interested in 
the youngest observed populations, we constrain our selection 
of stars to those with B-R < 0.5, which conespond to ages 
roughly ^ 500 Myr, acc ording to the Padova evolutionary 
models dGirardi et alj|2002l) . 

Fig.[T](middle) shows the spatial distribution of the selected 
sample of blue stars with B-R < 0.5. Since the final photo- 
metric catalog of dBW06 is the combined shallow photomet- 
ric catalog in the inner part with the deep catalog in the outer 
part of the galaxy, these authors checked the outer field for 
saturated stars that should have been incorporated in an equiv- 
alent shallow catalog of the outer field. They found that there 
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FIG. 2. — B-R, R CMD of all 208 894 stars detected with good photometr y in the observed S UBARU fields of NGC 6822. Indicative isochrones for 20 (blue) 
200 ( orange) and 500 Myr (red) fro m the Padova evolut ionary grid of models I Girardi et al. 2002) a re overlaid for a nominal metallicity Z = 0.004 ISkillman et al.l 
1989). reddening E(B-V) = 0.25 (Massey et al. 2007) and distance modulus m-M = 23.53 iGallart et al. 1996a). The green vertical line sets the limits of the 
"blue plume" stars, and the light blue box signifies the CMD area of the MS stars with (B-R) < 0.5, selected in this study for the identification of young stellar 
structures with ages down to roughly 500 Myr (see text in §[5J. [Color version of the figure will be available in the published electronic version of the paper.] 



are very few saturated stars present, and therefore the stellar 
distribution shown in the middle panel of Fig. Q] represents 
very well the whole observed blue stellar sample, despite the 
lack of a shallow photometry for the outer field. The cor- 
responding surface density distribution of the blue MS stars 
constructed by simple star counts in a smoothed grid of el- 
ements 64" x 64" in size is shown in the bottom panel of 
Fig.ffl 

From the maps of Fig. Q] one may conclude that the young 
blue stars (i) are distributed across the whole Hi disk, apart 
maybe from its southeastern part, (ii) they are mostly concen- 
trated in the central part of the galaxy in a "S-like" feature, 
and (iii) their distribution is rather clumpy. It should be noted 
that the density contour map of Fig. Q] (bottom) is a smoothed 
visualization of the large-scale stellar structures that exist in 
NGC 6822. The application, later, of the NN method (§© al- 
lows a far more detailed identification of the smallest possible 
stellar concentrations of the galaxy. 

Interestingly, the star counts of Fig. Q] revealed as an inde- 
pendent stellar structure the so-called "Northwestern Cloud", 
located at a position (R.A., Decl.) m (295.9,-14.6), which is 
speculated to be a separate system that is currently interact- 
ing with the main body of NGC 6822 and is thought to be 
responsible for the tidal arms in the southeast of the galaxy 
dde Blok & Walterll200d 120031) . It could also have triggered 
the star formation that eventually led to the creation of the 
large hole in the southeastern part of the main Hi disk, cen- 
tered at a position (R.A., Decl.) w (296.35,-14.9) (dBW06). 

The stellar content of the Cloud is dominated by a well- 
defined blue MS, and there are two low -brightness Hll regions 
identified in the Cloud (dBW06). Both findings indicate that 
the Cloud has experienced recent star formation. Unfortu- 



nately, the shallow SUBARU images do not cover its area, and 
therefore only stars up to mg ~ 20.5 are measured, the bright- 
est and bluest of which being well described by a logr ~ 7 
Padova isochrone. The presence, however, of a few saturated 
stars clustered in the same way as the unsaturated ones and 
the Ha observations suggest that even younger stars may be 
present in the Northwestern Cloud (dBW06). 

4. DETECTION OF CONCENTRATIONS OF YOUNG STARS 

The blue MS population selected in the previous section for 
the study of clustered star formation within the last ~ 500 Myr 
comprises 13 727 stars. In this section we apply the nearest 
neighbor (NN) method to investigate how these stars are clus- 
tered. A short introduct i on of the method, which i s des cribed 
by, e.g., ISchmeia et al.l ([2009), is given below (§ 14.11 ). Our 
method enables the detection of dense stellar concentrations, 
re veali ng the most prominent cluster population of NGC 6822 
(§ |4.2| i. While these systems span an order of magnitude in 
physical dimensions, we refer to all of them as clusters, al- 
though this term in principle applies only to a specific small- 
scale single-age class of systems, and not to the whole dis- 
tribution of detected objects. Independent repetitions of the 
method, for different density thresholds (§ [8]), revealed less 
dense stellar concentrations systematically belonging to larger 
ones, providing evidence of hierarchy in the distribution of the 
blue MS stars in NGC 6822. We refer to all these concentra- 
tions generally as structures. 

4. 1 . The Nearest Neighbor Density Method 

Star clusters are usually identified as regions of a cer- 
tain overdensity with respect to the background stellar den- 
sity. The NN method, introduced by Caserta no & Hull ( 1 19851) 
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FIG. 3. — The stellar density map of NGC 6822 (in logarithmic scale), constructed with the NN density method applied to stars with ages ^ 500 Myr for the 
10th NN. All stellar concentrations detected with density higher than 3<r above the average background are identified as members of the cluster population of 
NGC 6822. [Color version of the figure will be available in the published electronic version of the paper.] 



based on earlier work by Ivon Hoerneri d 19631) . estimates the 
local source density pj by measuring the distance from each 
object to its /'th nearest neighbor: 



1 



Pj : 



Sin) 



(1) 



where r, is the distance of a star to its jth nearest neighbor, 
S{rj) the surface area with the radius r, and m the average 
mass of the sources (m = 1 when considering number densi- 
ties). The chosen value of j depends on the sample size, and is 
correlated with the sensitivity to the density fluctuations being 
mapped. A small j value increases the locality of the density 
measurements while increasing the sensitivity to random den- 
sity fluctuations. On the other hand a large j value will reduce 
that sensitivity at the cost of some locality. 

The positions of the cluster centers are defined as 
the de nsity-weighted enhancement centers ( Case rtano & Hutl 
1985): 



(2) 



where xt is the position vector of the /th cluster member and p'j 
the jth NN density around this object. Similarly, the density 
radius rd en s is defined as the density-weighted average of the 
distance of each star from the density center: 



T,i \ x i- X d,i\p'j 
ZiP) 



(3) 



This ra dius typically correspo nds to the core radius of the 
cluster dCasertano & Hutll 19851) . 

We applied the NN method to our data considering severa l 
values for the nearest neighbors (j). ICasertano & Hutl (1 1985b 
have shown that low j values, in particular j = 1 or 2, are 
extremely sensitive to statistical fluctuations, therefore they 
suggest using a value of j > 6. On the other hand, the choice 
of a too large j value results in a loss of sensitivity to real 
density variations on smaller scales. Monte Carlo simulations 
have shown that a value of j = 20 is adequate to the detection 
of clusters with about 10 to 1500 members (B. Ferreira, pri- 
vate communication). However, since we aim at the detection 
of clusters even poorer than 10 members, we applied several 
test runs of the method that showed that j = 10 is a reasonable 
choice of number of NN. 

The 10th NN density map of the blue MS stars in NGC 6822 
is shown in Fig. [3] This map demonstrates that the selected 
blue MS stars are distributed in almost the whole extent of 
the disk of the galaxy and the northern-western cloud, as 
it has been demonstrated by the iso-density contour map of 
Fig. Q] The NN density map, however, along with the gen- 
eral (average) stellar distribution of these stars throughout the 
galaxy, can also track in detail individual stellar concentra- 
tions. Specifically, while the low-density distribution of the 
blue MS stars, shown in blue and green scales, clearly out- 
lines the disk of NGC 6822 and the north-western cloud, the 
high-density levels (shown in yellow and red) show a more 
clumpy and centrally concentrated stellar distribution, reveal- 
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ing individual stellar systems as density peaks. 

4.2. Detected Young Clusters 

A reasonable density level, above which the detected den- 
sity enhancements are accepted as star clusters is 3er above 
the average density, where a is the standard deviation of the 
background density. As a consequence, the original list of de- 
tected objects comprises all stellar density enhancements with 
density peaks at and above the 3a threshold. However, since 
our method does not set any dimensional limit for the detected 
objects, this list includes the smallest density peaks that could 
be possibly identified. Indeed, the smallest objects in the list 
are found to be minute density peaks, which include only one 
or two stars, and therefore they are most probably spurious 
detections at the level of the noise. Naturally, these detec- 
tions cannot be related to real physical stellar concentrations, 
but rather to the background density fluctuations. As a con- 
sequence, we exclude them from the final catalog of detected 
young clusters in NGC 6822, which thus comprises all stellar 
concentrations of 3 or more stellar members, that are found 
with our method to have 10th NN density values 3a above the 
average density. 

The detected clusters contain MS stars in the whole ob- 
served magnitude range of B 25. However, the range of 
brightness of the majority of the stars included in every clus- 
ter can be used to assign an indicative maximum age to each 
of them. We divided, thus, the observed catalog of MS stars 
into four magnitude ranges and we again computed the corre- 
sponding 10th NN density and applied a 3a threshold for the 
detection of clusters. Each of the selected magnitude ranges 
corresponds to a specific range of stellar ages, with the older 
age, specified by the turn-off of the best-fitting isochrone, be- 
ing associated to the faint magnitude limit. The four mag- 
nitude ranges were selected to include almost equal number 
of stars so that each run of the method would have the same 
statistical significance. We assign this age-limit to each of 
the ste llar sub-catalogs wit h the use of the evolutionary mod- 
els of iGirardi et al.l d2002l) for the distance, metallicity and 
reddening of NGC 6822, by applying isochrone fitting of the 
CMD. In this manner we are able to have a direct estimate of 
the age-range of the clusters identified in the corresponding 
stellar sub-catalogs, with those detected in the whole catalog 
having the oldest observable age of ~ 500 Myr. 

Although we cannot have an estimate of the actual age of 
each cluster, this approach is important for two reasons. First, 
age determination of the detected clusters by isochrone fit- 
ting of their individual CMDs is not really possible, due to 
their poor stellar number statistics. Moreover, the deriva- 
tion of ages from integrated surface brightness of each cluster 
would require the application of population synthesis tech- 
niques, which would introduce important model-dependent 
uncertainties, again due to low stellar numbers, especially 
for the faint clusters. The selected magnitude ranges are 
B < 21, B < 22, B < 23 and B > 23, which according to 
the models correspond to maximum ages for the clusters of 
Tmax — 95, 150, 250 and 500 Myr respectively. The majority 
of the detected clusters (36) are found to comprise stars in the 
faintest ranges or in the whole magnitude range and therefore 
r max ~ 500 Myr is assigned to them. Seven clusters are found 
with r max ~ 250 Myr, one with r max ~ 150 Myr, and three 
with r max ~ 95 Myr. The maximum ages of the clusters are 
given in column 4 of Table [3] for comparison to their dynam- 
ical timescales derived in § [7] In the following Section we 
discuss in detail the characteristics of the detected clusters. 
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FIG. 4. — Histogram of the size distribution of detected young star clusters 
in NGC 6822, with a bin size of 10 pc. Parameters of the best-fit Gaussian 
(red line) are also given. [Color version of the figure will be available in the 
published electronic version of the paper.] 

5 . THE CLUSTER POPULATION OF NGC 6822 

The catalog of identified star clusters in NGC 6822 consists 
of 47 clusters, the physical dimensions of which are defined 
by the 3cr density level in the 10th NN density map. Table Q] 
lists the identified clusters with the positions of their 10th 
NN density centers (Columns 2 and 3), the number of clus- 
ter members (Column 4), and their density and limiting radii 
(Columns 5 and 6 respectively) in seconds of an arc. The den- 
sity radius, rd ens , (Column 5) is derived from the NN method 
and Eq. (0 and roughly corresponds to the co re radius of each 
system. The so-called eff ective radius (e.g..lCarpenterJl2000h 
or equivalent radius (e.g.. iRoman-Zuniga et alJ l2008) of the 
systems (Columns 6 and 7) is the radius of a circle with the 
same area as the area enclosed by the cluster-defining 3a con- 
tour (r equ ; v = \/A c \/n). The total expected number of stars, N, 
(Column 10) and mass, M c \, (Column 11) are estimated from 
the extrapolation of the mass spectrum of each cluster assum- 
ing that it reassembles the average Galactic IMF (see § 15.21 ). 

5.1. Size Distribution of the Clusters 

From Table Q] it can be seen that our catalog of young clus- 
ters covers a variety of systems, starting with those of the min- 
imum of 3 stellar members with dimensions ^ 20 pc (cluster 
47), up to those with maximum size of almost 240 pc, includ- 
ing 77 stars (cluster 1). The size distribution of the clusters 
constructed by binning them according to their dimensions is 
shown in Fig. [4] Dimensions are given in physical units (pc) 
assuming a distance modulus of m—M = 23.53. A functional 
fit to this histogram shows that the size distribution of the de- 
tected clusters can be very well approximated by a normal dis- 
tribution. According to the best-fit Gaussian, drawn with a red 
line in Fig. [4] the dimensions of the detected systems are clus- 
tered around an average of ~ 68 pc with a standard deviation 
of ~ 18 pc. This is a rather interesting result, c onsidering that, 
as po i nted out earlier by various authors (e .g. jEfremov et alj 
1987; Ivanov 1996; Gouliermis et al. 2003), young stellar as- 
sociations in different galaxies seem to have typical dimen- 
sions close to this value (between 65 and 93 pc) with an av- 
erage of 80 pc. If, however, this length-scale is characteristic 
for star formation is stil l under debate (see,e.g., Basti an et al.1 
2007; Gouliermis 2010). 
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TABLE 1 

The Cluster Population of NGC 6822 with ages < 500 Myr, detected with 10th NN density 

VALUES 3(7 ABOVE BACKGROUND . 



Cluster 


R.A. 


Decl. 


«* 


'dens 


r, 


:quiv 




"'max 


N 


M a 


ID 


(deg) 


(deg) 




(") 


(") 


(pc) 


(mag) 


(M Q ) 




(10 3 M ) 


1 


a r\ r *~\ A f OO 

296.24692 


-14.74281 


77 


14.70 


50.29 


1 19.47 


19.34 


17.42 


87812 


O AO 1 O fQ 

9.99 ±0.68 


2 


296.22754 


1 A a oo 

-14.84399 


47 


10.24 


37.68 


89.50 


18.94 


19.33 


62032 


*1 1 f\ AC\ 

7.06 ±0.49 


3 


296.26373 


-14.91980 


29 


7.60 


29.18 


69.32 


20.35 


13.23 


47708 


r -1 a 1 r\ a 

5.43 ±0.38 


4 


296.25827 


-14.87879 


24 


7.09 


29.02 


68.93 


19.37 


17.27 


48364 


C CA 1 A OO 

5.50 ± 0.39 


5 


an/ anAllfl 

296.29990 


—14.88262 


24 


4.42 


25.76 


61.20 


21.19 


10.36 


29269 


OOO 1 A OO 

3.33 ± 0.23 


6 




-14.92559 


19 


5.50 


24.38 


57.91 


22.38 


7.17 


34399 


O A 1 I A O O 

3.91 ± 0.28 


7 


afi/; OOO 1 1 

296.25/1 1 


1 a onn 1 a 

— 14.89010 


15 


7.31 


21.20 


50.37 


21.89 


8.37 


12538 


1 HO I A 1 A 

1.43 ± 0.10 


8 


an/ a i i a 

296.26413 


1 a mn/ 1 

-14.93961 


15 


3.74 


20.20 


47.99 


21.52 


9.37 


34797 


O A/ 1 A OO 

3.96 ±0.29 


9 


OO/*"" aa 1 o/~" 

296.22186 


1 /I O 1 O 

-14.77319 


14 


4.04 


20.10 


47.76 


20.19 


13.82 


19757 


O O C 1 A 1^ 

2.25 ±0.16 


10 


•~w\/ oaoo 1 

296.20294 


1/1 'O O O /I /I 

-14.73244 


16 


3.84 


19.99 


47.48 


19.30 


17.58 


16055 


1 O O I A "1 O 

1.83 ±0.13 


1 1 


a n 00^7 1 a 

296.23712 


1 /i Taoin 

-14.72879 


14 


4.39 


18.20 


43.23 


20.25 


13.61 


16998 


1 AO I A 1 A 

1.93 ±0.14 


12 


O O /~ O O O C /I 

296.29254 


1 A O O "7 O A 

-14.88720 


1 1 


4.26 


17.64 


41.90 


21.92 


8.29 


18219 


O /Y~7 1 A 1 C 

2.07 ± 0.15 


13 


an/ ioio 1 

296.28/81 


1 A OO £LT~i 

— 14.936/3 


10 


6.29 


17.3 1 


41.12 


20.03 


14.46 


26106 


O A "7 1 A OO 

2.9/ ± 0.22 


14 


an/ TJ /T 1 1 

296.2361 1 


1/1 O A A C A 

—14.84454 


10 


2.81 


17.24 


40.94 


21.49 


9.46 


21508 


O A C 1 A 1 O 

2.45 ±0.18 


15 


anil Ti oao 

296.31003 


— 14.92668 


12 


2.13 


17.20 


40.86 


22.16 


7.68 


25761 


O AO 1 A OO 

2.93 ± 0.22 


16 


296.24265 


1 a "7 con A 

— 14.75290 


10 


2.99 


16.30 


38.73 


20.82 


1 1.55 


21005 


O OO 1 A 1 "7 

2.39 ± 0.1 / 


17 


aril' ai-ia i a 

296.29214 


1 /i 00 "7 1 

— 14.90718 


13 


5.14 


16. 17 


38.40 


21.69 


8.89 


10096 


1 1 C I A OO 

1.15 ± 0.08 


18 


a r\ s a a f \ r a 

296.22964 


-14.73439 


10 


6.30 


16.01 


38.03 


20.95 


11.14 


10703 


1 00 1 a 00 

1.22 ± 0.08 


19 


296.22729 


1/1 T1 Oil 

-14.71831 


1 1 


1.67 


16.01 


38.03 


21.43 


9.64 


22769 


O CO 1 A "1 O 

2.59±0.19 


20 


ari/ anaaf 

296.2/225 


1/1 OO COO 

—14.88502 


10 


3.15 


15.72 


37.34 


21.09 


10.67 


34637 


O A A 1 A OA 

3.94 ± 0.30 


21 


296.21454 


1 A "7*£ AO /i 

— 14.76024 


10 


4.29 


15.47 


36.75 


21.23 


10.25 


22983 


O H 1 I A 1 O 

2.61 zb 0.19 


22 


^\C\f a a o a r 

296.23825 


1 /I 'O O ^ O A 

-14.73670 


9 


4.81 


15.46 


36.73 


20.75 


1 1.80 


12147 


1 O O I A 1 A 

1.38 ± 0.10 


23 


aril' O iCTA*! 

296.26/06 


1 a rnnrin 

— 14.92909 


7 


3.51 


15.27 


36.28 


20.81 


1 1.60 


20109 


O OO I A 1 "7 

2.29 ± 0.17 


Z4 


296.26959 


1 /i 0*7000 
— 14.8/903 


9 


4.01 


14.69 


34.89 


22.70 


0.4j 


14426 


1 1 A 1 O 

1.64 ± 0. 12 


25 


296.21567 


1/1 O /I A CC 

— 14.84455 


8 


4.55 


14.26 


33.88 


21.36 


9.85 


31207 


O CC 1 A 0"7 

3.55 ± 0.27 


26 


an/' 1 O 1 O O 

296.13132 


1 /i /Afire 

-14.69955 


1 1 


2.95 


14.06 


33.39 


19.58 


16.35 


38558 


/IOO 1 A O C 

4.39 ±0.35 


27 


296.26364 


-14.90800 


8 


2.39 


13.98 


33.22 


21.55 


9.30 


29133 


3.31 ±0.26 


28 


aril' OOCO/1 

296.2/594 


1 /I OO /I OA 

— 14.92400 


7 


4.27 


13.55 


32.18 


21.59 


9.18 


9849 


1 1 O 1 A OO 

1.12 zb 0.08 


29 


296.22586 


— 14.7 /637 


7 


3.80 


13.03 


30.95 


20.43 


12.94 


25421 


O OO 1 A OO 

2.89 ± 0.22 


30 


296.2 / 1 18 


1 A on 1 QO 

— 14.89139 


7 


3.03 


12.82 


30.46 


23.03 


5.79 


17718 


i n 1 _i_ a 1 c 
2.01 ± 0. 15 


3 1 


an /- <~s£L A on 

296.26489 


1/1 ot 1 / a 

—14.87163 


4 


1.76 


12.80 


30.42 


20.93 


1 1.20 


1 1498 


10 1 I A OO 

1.31 ± 0.09 


32 


296.22665 


1 /I OT7T7 

— 14.87/27 


9 


2.60 


12.79 


30.38 


21.38 


9.79 


26342 


O OA 1 A OO 

3.00 ± 0.23 


33 


296.22006 


1/1 /i ic ~7 n 

— 14.846/2 


6 


2.38 


12.67 


30.09 


22.24 


7.50 


31063 


O CO 1 A OO 

3.53 ± 0.28 


34 


a n /~ OO^T 1 \ 

296.2271 1 


1 A aaa 1 r 

-14.72315 


8 


2.98 


12.33 


29.30 


20.86 


1 1.41 


26343 


O O A 1 A O O 

3.00 ±0.23 


35 


296.20230 


-14.76187 


6 


4.90 


11.71 


27.81 


20.88 


11.37 


13942 


1.59±0.12 


36 


296.19022 


-14.85646 


7 


2.68 


11.49 


27.29 


19.84 


15.22 


18638 


2.12±0.16 


37 


296.31396 


-14.97089 


5 


2.49 


11.33 


26.90 


21.46 


9.55 


15824 


1.80±0.14 


38 


296.20490 


-14.75966 


5 


2.06 


11.11 


26.40 


21.70 


8.88 


9558 


1.09 ±0.08 


39 


296.23492 


-14.72453 


6 


2.86 


10.97 


26.06 


20.87 


11.38 


12427 


1.41 ±0.11 


40 


296.23090 


-14.75460 


6 


2.81 


10.37 


24.64 


20.92 


11.22 


17777 


2.02 ±0.16 


41 


296.23541 


-14.73430 


4 


2.33 


9.63 


22.87 


22.36 


7.21 


10157 


1.15 ±0.09 


42 


296.19272 


-14.86216 


3 


2.43 


9.00 


21.39 


22.82 


6.22 


12311 


1.40±0.11 


43 


296.24420 


-14.73231 


4 


1.62 


8.27 


19.65 


22.86 


6.13 


13818 


1.57 ±0.12 


44 


296.31229 


-14.75217 


4 


1.99 


7.41 


17.59 


21.29 


10.05 


11492 


1.31 ±0.10 


45 


296.25198 


-14.88950 


4 


2.42 


5.78 


13.74 


22.49 


6.91 


20460 


2.32 ±0.20 


46 


296.24365 


-14.87206 


3 


2.52 


5.58 


13.25 


21.42 


9.66 


14739 


1.68±0.14 


47 


296.21973 


-14.88315 


3 


4.02 


3.85 


9.14 


22.18 


7.65 


8670 


0.99 ±0.09 



NOTE. — A detailed description of the parameters is given in §[5] 



5.2. Masses of the Clusters 

The sample of detected clusters naturally includes a variety 
of systems not only concerning their physical dimensions but 
also their stellar content and in consequence their total mass. 
While our data are not deep enough to cover the whole extent 
of masses for the stellar members of the detected clusters, and 
not complete enough to provide us with the actual numbers 
of stars (per luminosity or mass range) for each of them, it 
is worthwhile to perform a first-order calculation of the total 
stellar mass included in each detected cluster, based on sev- 
eral assumptions. We identify the number of stars included in 
each cluster and their luminosities from our photometric cata- 
log. Then, in order to have a first estimation of the total mass 
of each cluster we (i) build the observed Luminosity Func- 
tion (LF) of each cluster, (ii) apply a mass-luminosity relation 
(M-LR) based on the stellar evolutionary models for the con- 
version of this LF to the Mass Function (MF) of the cluster, 



and (iii) extrapolate this observed MF to the lower masses for 
the calculation of the total mass of each cluster. 

5.2.1. Mass-Luminosity Relation 

For the establishment of the conversion of stellar lu- 
minosities to masse s, we use the evolutionary models by 
iGirardi et alJ ((2002). Our clusters identification is applied 
for all detected stars of NGC 6822 formed within the last 
~ 500 Myr, and therefore the identified clusters are expected 
to have different ages from each other within this time range. 
As a consequence, in order to construct a realistic M-LR for 
each cluster, one should establish a correct age with the use 
of the CMD of the detected stars in the cluster. However, the 
construction of individual M-LR for each cluster is not pos- 
sible, because in all cases the stellar numbers found within 
the detected clusters are not sufficient to build complete and 
informative CMDs. 
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FIG. 5. — Mass-Luminosity relation c onstructed w ith the use of isochrones 
from the grid of evolutionary models by Girardi et al. 1 2002) for the metallic- 
ity and at the distance of NGC 6822. Isochrones of a large variety in ages are 
used for the correct quantification of the uncertainties in the determination of 
the cluster masses due to age differences among the identified clusters. [Color 
version of the figure will be available in the published electronic version of 
the paper.] 

On the other hand, one may construct a global M-LR to be 
used for the conversion of LFs to MFs for all clusters, but this 
M-LR should be established from evolutionary models that 
cover a large range in ages. Indeed, as it is seen in the CMD of 
Fig.|2]the observed MS is populated by stars in an increasingly 
larger variety of ages towards fainter magnitudes down to B ~ 
25 mag. Therefore, we construct a global M-L R for all iden- 
tified c lusters using the evolutionary models of Girar di et alj 
(2002) not for specific ages, but for all available ages and for 
the appropriate metallicity of Z = 0.004. With this treatment 
the uncertainties in stellar masses due to the age differences 
among the identified clusters will be inherited to the global 
M-LR and affect accordingly the determination of the total 
mass of the clusters. The derive d M-LR, for an assu med av- 
erage reddening E(B-V) = 0.25 ( Massev et aTll2007h and dis- 
tance modulus m—M = 23.53 (Ga llart et al.lll996ah . is shown 
in Fig. [5] From this figure it can be seen that indeed the stellar 
mass uncertainty is a function of the MS magnitude, with the 
larger uncertainties appearing for the parts of the MS where 
models of different ages show the greatest differences in the 
derived masses. While B ~ 25 mag corresponds more or less 
to our detection limit, and our stellar sample does not include 
stars brighter than B ~ 17.5 mag, we construct our M-LR for a 
more extended magnitude range for reasons of completeness 
in its shape. This M-LR is being used for the direct trans- 
formation of the luminosity of each star detected within the 
limits of one of the NGC 6822 clusters into mass. 

5.2.2. The Stellar Mass Spectrum of the Clusters 

The stellar mass function is the number distribution of stars 
according to their masses, constructed for a given volume of 
space in a cluster. This function is known as the stellar Present 
Day Mass Function, which is usually called simply the Mass 



Function of the cluster. In this study we refer to this function 
as the Stellar Mass Function (SMF), to distinguish it from the 
Cluster Mass Function discussed later. The SMF of a cluster is 
described by the function £(logm), which gives the number of 
stars per unit logarithmic (base ten) mass interval d login per 
unit area (e.g. /kpc 2 ). Alternatively, one may refer to the Stel- 
lar Mass Spectrum (SMS) of the cluster, f(m), which is the 
number of stars per unit mass interval dm per unit area. The 
common use of both these functions is based on the defini- 
tion of the stellar Initial Mass F unction (IMF), vari ous forms 
of whi ch are discussed by, e.g., Kroupa (2002) and Chabrier 
(2003). All these distributions are usually parametrized by 
their indices. Specifically, for the SMS f(m), and the SMF 
£(logra), these indices are defined as 

<flog/(m) 

7 = — 7j (4) 

d log m 

and 

r ^ aflogcl(logffl) 
dlogm 

These two indices are basically the logarithmic derivatives or 
slopes of f(m) and £(logm) respectively an d for power -law 
distributions they are independent of mass (Scald 1 19861) . A 
reference value for the SMF slope is T = -1.35, which is the 
index of the classical IMF for stars in the solar neighborhood 
with masses 0.4 < M/M Q < 10, found bv lSalpeterl (11955b . 
The corresponding SMS index is j = T— 1 ~ -2.35. For 
comparison, the lognormal field-star IMF by Mill er & Scalo I 
( 119791) has T ~ -(1 +logm). A basic relation between SMF 
and SMS i s £(logm) = (In 10) • m ■ f(m) ~ 2.3 • mf(m) (see 
Scalo Fi986l) . 

The construction of the SMS of each cluster is essential for 
the calculation of its expected total mass through extrapola- 
tion of its SMS. However, in most of the cases of identified 
clusters in NGC 6822 the detected stellar numbers are not 
sufficient for the construction of a meaningful SMS for each 
cluster. This problem can be surpassed with the use of a com- 
mon SMS constructed by all stars detected within the identi- 
fied clusters. This approach, naturally, assumes that the SMS 
is indeed univer sal, in line with most findings in th e local uni- 
verse (see, e.g., Massey 2006; Bastian et al. 2010). We con- 
struct the global SMS of the clusters in NGC 6822 by count- 
ing all stars detected within the limits of all the 47 identified 
clusters according to their masses, and distributing the stellar 
masses in mass-bins 1 M wide. This SMS is given in Ta- 
ble |2j and it is shown in Fig. [6] This figure also shows the 
best-fit line derived from the application of a linear regres- 
sion to the most complete mass-bins corresponding to masses 
between ~ 5 and 13 M . 

The determination of the slope of power laws from uni- 
formly binned data using linear regression is found to 
comprise biases, which are caused by the correlation be- 
tween the number of stars per bin and the assigned weights 
dMafz Apellaniz & Ubeda 2005). However, we measured the 
average mass of the stars counted in each bin for the con- 
struction of the global SMS, and we found that it is almost 
equal to the mean mass of the bin. The corresponding stan- 
dard deviations per bin, are found to be smaller than 1 M Q , 
the width of each bin, and therefore there is no "leak" of stars 
expected from one mass-bin to its neighboring, and thus no 
random fluctuations in the constructed SMS. The measured 
average mass per bin, rather than the mean mass of each bin, 
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TABLE 2 
THE GLOBAL Stellar Mass Spectrum 
CONSTRUCTED FROM ALL STARS DETECTED 



WITHIN THE LIMITS OF THE 47 IDENTIFIED 




CLUSTERS IN NGC 6822. 




Altars 


(m) 




P* 


Mq 




M Q 


M Q 


stars/kpc 2 /M0 


2.50 


1 


2.94 




7.43 ±3.71 


3.50 


31 


3.75 


0.03 


118.83 ±20.68 


4.50 


123 


4.51 


0.02 


460.48±41.19 


5.50 


124 


5.46 


0.04 


464.19±41.35 


6.50 


84 


6.42 


0.01 


315.65 ±34.04 


7.50 


59 


7.35 


0.04 


222.81 ±28.52 


8.50 


45 


8.49 


0.05 


170.82 ±24.91 


9.50 


31 


9.41 


0.04 


118.83 ±20.68 


10.50 


16 


10.50 


0.11 


63.13 ±14.85 


11.50 


15 


11.36 


0.04 


59.42 ±14.38 


12.50 


6 


12.52 


0.17 


25.99±9.10 



NOTE. — The mean mass per bin, given in Col. 1, 
is comparable to the measured average mass of Col. 
3, calculated from all stars counted per bin, given in 
Col. 2. The corresponding standard deviations, given 
in Col. 4, do not exceed the width of each bin (1 Mq). 
The counted stellar numbers per surface used for the 
construction of the SMS, and the corresponding Pois- 
son errors are given in Col. 5. Values are plotted in 
Fig.[6] 
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FIG. 6. — The Stellar Mass Spectrum of all stars comprised within the areas 
covered of all 47 detected clusters with age $ 500 Myr in NGC 6822. The 
total sample includes 561 stars with masses between ~ 3 and 13 Mq. The 
sample is complete for stars down to ~ 5 - 6 Mq . A linear fit (solid line) to 
the SMS down to this mass limit gives a slope of 7 between -2.21 and —2.82, 
very close to the typical Galactic field slope for the same mass range. The 
plotted points correspond to the average mass of all stars counted per mass- 
bin with the x-axis errors derived from the standard deviation of the masses 
p er bin. The errors on the y-axis reflect the Poisson statistics. See also text in 

§11221 

is plotted in Fig. and the corresponding <7/ m \ are shown as 
error bars on the x-axis, demonstrating the accuracy of the 
pointing of each mass-bin. The linear regression for the most 
complete mass range of 5 ^ m/M Q ^ 13 provides a SMS 
slope 7, which is exceptionally close to the average Galactic 
field value, equal to -2.5 ± 0.5 (black-solid line in Fig. [6]). 

5.2.3. Total Masses of the Clusters 

An estimate of the total mass of each cluster in NGC 6822 
can be achieved through integration of its SMS. However, 
as we discuss earlier, the construction of a meaningful SMS 



for each cluster is not possible due to low stellar numbers, 
and therefore we make use of the global SMS assuming that 
this distribution is universal and representative of all clusters. 
While, according to our assumption, the slope of the SMS 
should not be different from one cluster to the other, the abso- 
lute stellar numbers per mass-bin of the SMS in each cluster 
are different and thus the corresponding total cluster mass. In- 
deed, the individual SMS of the clusters do deviate from the 
global SMS, mostly in stellar numbers and, in some cases of 
clusters with very few stars, in the slope. This is demonstrated 
in Fig.|7l where the SMS of seven selected indicative clusters 
from our sample are plotted. In this figure it can be seen that 
the general trend of the SMS slope indeed seems consistent 
with the global SMS (solid line), but the stellar numbers (per 
surface unit) are quite different. These differences will deter- 
mine the differences in total mass of the clusters. We estimate 
the cluster masses in three steps: 

(i) We first extend the observed global SMS to the low- 
mass regime below our detection limit according to the 
generally accepted Galactic field mass spectrum. This 
mass spectrum is usually parametrized with a series of 
power laws, with exponents c hanging in different mass 
ranges (see, e.g., IScaloll 19981). The most recent param- 
eterization is that by Kroupa ( 2002), according to which 
the slope 7 changes from 7 = -0.3 in the substellar mass 
range, to 7 = —1.3 for masses between 0.08 M Q and 
0.5 M , 7 = -2.3 for 0.5 < m/M© < 1, and 7 between 
-2.7 and -2.3 (±0.3) for stars of higher masses. This 
SMS is generally characterized as the Galactic average, 
in the sense that it is reasonably valid for different re- 
gions of the Galaxy. 

Considering that the SMS of all clusters identified in 
NGC 6822 should have a multi-power law form, with 
exponents that change at the same mass limits as the 
Galactic SMS, we assume that the global SMS of the 
cl usters has the average slope of those measured in 
§ 15.2.21 equal to 7 = -2.5 ± 0.5, for stars with masses 
down to our completeness limit of about 5.5 M and 
we extrapolate it with the same slope down to the mass 
limit of 1 M . For the extrapolation of the global SMS 
to sub-solar masses we consider the exponents of the 
Galactic SMS, as they are discussed above for stars with 
masses down to ~ 0.1 M . 

The extra polated global SMS, following the parameter- 
ization by Kroupa ( 2002), is thus a three-part power law 
function of the form 

/(/«)« ^— with/= 1,2,3. (6) 

The slopes 7 depending the mass range are 

71 =-1.3 ±0.5 , 0.1 <m/M <O.5 

72 = -2.3 ±0.3 , 0.5 <m/M Q <1.0 (7) 

73 =-2.5 ±0.5 , 1.0 <m/M Q . 

(ii) We then normalize the extrapolated global SMS to fit 
the stellar numbers observed in every cluster. Specifi- 
cally, we normalize the extrapolated global SMS to the 
mass-bin of around 5 to 6 M Q , because all detected 
clusters, including those with only three stars, show to 
have at least one star at this specific mass range. There- 
fore, by normalizing the global SMS to this mass limit 
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FIG. 7. — Indicative Stellar Mass Spectra (SMS) from the sample of clus- 
ters detected in NGC 6822. Different symbols are used for plotting the SMS 
of seven indicative clusters, comprising a variety of objects from the richest 
cluster in the sample (Cluster # 1 ) plotted with filled circles to one of the poor- 
est (Cluster # 41) with only two mass-bins plotted with open boxes. Points 
for individual clusters are also connected together with lines. The slope of the 
global SMS is represented by the solid line, as in Fig. [6] This graph demon- 
strates the variability in stellar numbers among the clusters, and the inability 
to extrapolate the observed SMS of each cluster, due to low stellar numbers 
and thus the necessity for the use of a global SMS. The mass range over which 
the global SMS is considered for extrapolation for all clus ters, m ~5-6Mq, 
is indicated by the vertical shaded area. See also text in § 15.2.31 

we apply a comparable approach to all, rich and poor, 
clusters. 

(iii) We finally integrate the normalized SMS of each clus- 
ter from the largest mass observed in the cluster to the 
smaller stellar mass assumed to exist in all clusters of 
0.1 Mq, and derive an estimation of the total cluster 
mass, from its surface stellar density 

M c i = p* 7T rg quiv , (8) 

where 

p*= / f(m)dm. (9) 
Jo.i 

The total stellar numbers of the clusters derived from this 
method and the corresponding total stellar masses are given 
in Table Q] (Cols. 10 and 11 respectively). 

6. MASS DISTRIBUTION OF THE CLUSTERS 

The total stellar masses of the clusters of NGC 6822, as de- 
rived above, cover a range over an order of a magnitude with 
the smaller cluster mass being around 10 3 M . The num- 
ber distribution of the clusters according to their masses, i.e. 
the mass spectrum of the clusters, which is widely termed the 
Cluster Mass Function (CMF) is shown in logarithmic scale in 
Fig. [8] The CMF derived from young ( ^ 10 Myr) star clus- 
ters is defined as the Cluster Initial Mass Function (CIMF). 
Considering that the CMF follows the functional form of a 
power-law, as it is generally assumed, we applied a linear re- 
gression to the data of Fig. [8] The best fit returns the function 

log/V cl = (-1.47 ±0.72) logM cl + (1.60 ±0.31). (10) 

We found, thus, that the CMF of the massive clusters of 
NGC 6822 can be described by a power law with index 
~-1.5±0.7. 



6.1. Comparison with the CIMF 

Previous works on the CIMF in a variety of galactic envi- 
ronments have shown that this func tion is well d escribed by a 
power law with index -2 (see, e.g., Gieles 2009): 

dN 

— «M-«, a = 2. (11) 

Specifically, the CMF of star clusters in the Antennae galaxies 
(NGC 4038/9) with ages ^ 160 Myr is found, over the range 
10 4 $ M cl /M w < 10 6 , to have an i ndex of -2 dZhang & Falll 
119991) . iMcCrady & Graham! (120071) constructed the CMF of 
19 super star clusters (SSCs) with M cl > 10 5 M Q and ages 
of the order of 10 Myr in the nuclear starburst of M 82, and 
found a power-law with an index -1.91 ±0.06. The SSCs 
with an average age ^10 Myr in the nearby starb urst galaxies 
NGC 3310 and NGC 6745 are investigated by Ide Griis et alj 
(2003), who find CMF indexes of -2.04 ± 0.23 and -1.96 ± 
0.19 respectively. These CMFs correspond to the CIMF due 
to the youthfulness of the cluster samples. All three studies on 
massive star clusters agree to an almost identical CIMF index 
of -2. In addition, a similar CIMF index of -2.1 ±0.3 is found 
for 10 Myr old clusters in the mass range 2.5 10 3 < M c i/M < 
5 10 4 in the inner spiral arms of M 5 1 dBik et alj2003l) . All the 
aforementioned studies provide evidence of a "universality" 
of the CIMF. Since most of the clusters in our sample have 
undergone substantial dynamical and stellar evolution, their 
CMF cannot be compared with the CIMF. It is interesting, 
though, to notice that the derived CMF index of NGC 6822 in 
the mass range 10 3 - 10 4 M Q of ~ -1.5 ±0.7 is well within lcr 
of the "global" CIMF, as the large uncertainties in our CMF 
do not allow us to identify any significant difference from it. 

6.2. Comparison with the Magellanic Clouds 

Our CMF of NGC 6822, which covers stellar concen- 
trations with ages of up to r ~ 500 Myr, should have 
been significantly affected by evolutionary effects and there- 
fore it is not directly comparable to the CIMF. Compara- 
ble CMFs in terms of cluster ma sses and ages to t h at de- 
rive d here are those pres e nted b y Ide Griis & Anders! ( 120061) 
and Ide Griis & Goodwin] d2008l) for the cluster populations 
in the Large and Small Magellanic Clouds (LMC, SMC) re- 
spectively. The CMF index in the LMC is given as a func- 
tion of minimum ma ss and maximum age, based on mass- 
limited samples (see Ide Griis & An ders 2006, their Table 3 
and Fig. 8). For log(m cl /M ) min = 3 and logfr/yr) =8.7 
which corresponds to our sample. Ide Griis & And ers (2006) 
find a CMF index -1.98 ±0.08, an average value which is 
somewhat different than ours, but nevertheless covered by our 
uncertainti es, and quite similar to t he CIMF. In addition, in 
the SMC Ide Griis & Goodwinl (12008b find a CMF with index 
a ^ 1.2 for almost the same cluster mass range with ours. 
This slope falls within lcr of our CMF slope, and is therefore 
also consistent with it. 

7. STRUCTURAL PARAMETERS OF THE CLUSTERS 

We estimate the stellar density, q c \, in the half-mass radius, 
n,, of each cluster and the disruption time, t&, of the cluster due 
to inter action with passing-by in terstellar clouds, as described 
in, e.g., [Gouliermis et al. (2002, their section 6). The former 
is calculated assuming spherical symmetry for all clusters as 

_ M cl /Mq \ 

Qa = 4 — 7 —5- , (12) 
§7rr 3 VPC 3 / 
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FIG. 8. — The Mass Spectrum of the star clusters or the Cluster Mass Func- 
tion of NGC 6822, based on calculation of their total stellar masses in § 15.2.31 
under the assumption that all clusters follow a common high-mass stellar 
IMF, extrapolated to the sub-stellar limit according to the Average Galactic 
IMF. [Color version of the figure will be available in the published electronic 
version of the paper.] 

while the latter is given as 



t d = 1.9 x 10 gel (years). 



(13) 



The dynamical status of a young star cluster is defined by 
two additional time-scales, the crossing and the two- body re- 
laxation time, which are given as (e.g., Kroupa 2008) 



2r h 

— and 

a 



frela 



N 

InN 



(14) 



respectively. The three-dimensional velocity dispersion of the 
stars in the cluster, a, is given as 



GMci 



(15) 



where e is the star formation efficiency (SFE) and the half- 
mass radius of the cluster. We estimate these time-scales for 
the clusters of our sample, with the application of Eqs. (fl4l . 
assuming the same SFE with several nearby Galactic gas- 
embedded clusters, which has been fou nd to range typically 
from 10% to 30% dLada & Lada 1 120031) . In our calculations 
we assume an average value of e ~ 0.2, which leads to values 
of a between 12 and 89 km s" 1 for our clusters. 

For the measurement of the rj, of each cluster we estimate 
the so-called Spitzer radius, rs p . Considering that ~ 0.9 rs p 
only for dynamically relaxed spherical clusters (e.g., Spitzer 
1969), we can only approximate the actual of our clusters. 
The Spitzer radius is a dynamically stable radius, defined by 
the mean-sq uare distance of the stars from the center of the 
cluster (e.g., Spitzer 1987|) as 



■ J i=i'i 
N 



(16) 



where r, is the projected radial distance of the /th star of the 
cluster in a total sample of N stars. We estimate this radius for 
each cluster in our sample only from the stars identified within 
its limits, for which there are positions available, and not from 
the expected total number of stars according to the extrapo- 
lated mass spectrum of each cluster. As a consequence, the 
accuracy of the derived rs p of each cluster is subject to the 
number of its identified stars. 



m 2.5 



J 2.0 




0.5 



3.5 



the cluster population of NGC 6822 (10th NN, 3a density detection level). 
This c orre lation is consistent with that found for Galactic embedded clusters 
(see §|7T}. 

The estimation of rs p allows the evaluation of u and g c \ of 
the clusters and consequently of their t&, t cr and f re i ax according 
to equations ( fT3l ) and (TBI) . For the latter we use the expected 
total number of stellar members, N, as it is derived from the 
extrapolation of the mass spectrum of each cluster. All addi- 
tional estimated structural parameters for the young clusters 
of NGC 6822 are provided also in Table [3] 

7.1. Parameters correlations of the clusters in NGC 6822 

A correlation between the number of stars in a cluster and 
the radiu s of the cluster has been found in Ga lactic embedded 
clusters dAdams et alJl2006tlAllen et alj|2007l) : While n* and 
r vary by about two orders of magnitude, the average surface 
density of cluster members (n^/nr 2 ) is nearly constant. We 
find the same correlation, following a power law, between n* 
and r equ iv for the clusters identified in NGC 6822, which span 
an even wider range in radius, n*, and age (Fig.|9|l. The corre- 
lation has a sl ope in the log-log plo t of k ~ 0.6, comparable to 
that found by Adam set al.l d2006l) for their sample of young 
Galactic clusters. Our fit, which has a correlation coefficient 
of q ps 0.92, is affected by the smallest clusters (n* < 5) that 
show some scatter. A similar correlation exists also between 
n* and fdens, but it becomes less clear when using the expected 
number of stars N instead of the observed number A cor- 
relation between number of stars and radius is also observed 
for the stellar structures we identify below with the applica- 
tion of the NN method with density levels different than 3cr 
(see §H3). 

8. HIERARCHICAL CLUSTERING IN NGC 6822 

The majority of stars fo rm in clusters and a g gregates of 
variou s sizes and masses (Lada & Lada 1 120031 : lAllen et alj 
2007). Young star clusters are often found inside larger 
complexes, which can be parts of even larger structures. In 
turn, stellar clusters often consist of distinct subclusters, ap- 
pearing to for m a hierarchy of systems over a wide range 
of sca les (e.g. lEfremov & Elmegreer] [T998: Elmegree n et alj 
2000). The ISM is also hierarchically structured (sometimes 
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TABLE 3 

Structural parameters of the star clusters of NGC 6822. 



Cluster 


R.A. 


Decl. 




Qd 




a 


fd 




^relax 


ID 


(deg) 


(deg) 


(Myr) 


(M /pc 3 ) 


(pc) 


(kms-') 


(Myr) 


(Myr) 


(Myr) 


1 


TA/' O /I /~"AO 

296.24692 


—14.74281 


500 


A AC 1 A AA 

0.05 ± 0.00 


36.92 


89.0 


9.00 


0.81 


626.83 


2 


296.22754 


1/1 o a o aa 

—14.84399 


500 


A AA 1 A A 1 

o.o9±o.oi 


26.87 


63.8 


16.49 


0.82 


463.66 


3 


^M~\s o/"o^O 

296.263 73 


1 a n i noA 

—14.91980 


500 


A 1 A I A A 1 

0.14±0.0l 


21.16 


49.7 


25.97 


0.83 


369.59 


4 


mi; o c o o o 

296.2582/ 


1 /I 07070 

—14.8 /8 /9 


500 


A 1 1 I A A 1 
0.1 1 ±0.01 


22.79 


51.9 


2 1 .07 


0.86 


385.82 


5 


296.29990 


1 /I o oozio 

—14.88262 


500 


A O O 1 A Al 

0.33 zb 0.02 


13.48 


31.1 


61.66 


0.85 


242.05 


6 


oa/~" oA/1ao 

296.29492 


—14.92559 


500 


0.21 ±0.02 


16.40 


37.1 


40.14 


0.87 


285.19 


7 


296.20/1 1 


1 /i o 1 a 

—14.89010 


500 


A AC 1 A AA 

0.05 zb 0.00 


18.43 


23.8 


10.32 


1.52 


201.97 


8 


296.26413 


—14.93961 


500 


A £L"i 1 A AC 

0.63 ± 0.05 


1 1.42 


31.2 


120.41 


0.72 


OOO AO 

238.9/ 


9 


296.22186 


1 A T7T 1 M 

—14.77319 


500 


A OA 1 A AT 

0.39 ±0.03 


11.11 


23.2 


74.29 


0.94 


187.72 


10 


anon 1 

296.20294 


1/1 T~}<~S A A 

—14. 73244 


150 


A OA 1 A AO 

0.30 ±0.02 


1 1.28 


21.0 


57.75 


1.05 


174.13 


1 1 


^ f"\ O O 1 1 O 

296.23712 


i /i ^toota 

— 14.72879 


250 


A It 1 A A>"l 

0.25 ± 0.02 


12.31 


22.6 


47.05 


1.07 


186.06 


12 


296.29254 


1 A O O 1 O A 

— 14.88720 


500 


A OO 1 A AO 

0.32 ± 0.02 


1 1.54 


22.7 


61.08 


1.00 


185.32 


13 


o a/: ooto 1 

296.28/51 


1 yl A O £ T O 

— 14.936/3 


500 


A 1 O I A A 1 

0.18 zb 0.01 


15.92 


31.9 


33.31 


0.98 


251.41 


14 


o A/; oo /; 1 i 

296.2361 1 


1/1 O A A C /I 

—14.84454 


250 


1 OC I A A A 

1.25 zb 0.09 


7.76 


20.2 


237.77 


0.75 


162.28 


15 


296.31003 


—14.92668 


500 


2.16 zb 0.16 


6.87 


20.8 


410.02 


0.65 


164.18 


16 


296.24265 


1 A OCOAA 

—14. 75290 


250 


A A"7 1 A A"7 

0.97 zb 0.07 


8.37 


20.7 


184.95 


0.79 


166.92 


17 


296.29214 


1 A Am 1 o 

—14.90718 


500 


A A"7 1 A A 1 

0.07 zb 0.01 


15.45 


19.5 


14.13 


1.55 


169.78 


18 


O A/C q/in/ a 

296.22964 


1 A TO A Tfl 

—14.73439 


95 


A C\C _l_ A AA 

0.06 zb 0.00 


16.53 


20.8 


12.23 


1.56 


179.65 


19 


296.22729 


1/1 H0O1 

— 14.71831 


500 


O A O 1 A o*r 

3.48 ±0.26 


5.62 


17.7 


660.40 


0.62 


141.35 


20 


296.2/225 


i/i oo cao 

—14.88502 


500 


1 A A I A 1 1 

1.44 zb 0.1 1 


8.67 


27.1 


OOO OO 

215.13 


0.63 


OAO OO 

20/. 82 


21 


O A/C ai j 

296.21454 


1 A O £AO A 

— 14. 76024 


500 


A O £L 1 A AO 

0.36 zb 0.03 


1 1.97 


25.9 


69.20 


0.90 


206.96 


22 


OA/; O O O O C 

296.23825 


1 A 1 1 /LH\ 

—14. 73670 


500 


A 1 A I A A 1 

0.14zb 0.01 


13.46 


20.0 


25.71 


1.32 


170.37 


23 


OA/; <~S£L~1(\£. 

296.26/06 


i /i nnnan 

—14.92909 


95 


A O 1 I A AC 

0.71 zb 0.05 


9.15 


21.2 


135.34 


0.85 


171.56 


24 


o a/: T/:ncn 

296.26959 


1 yl 0"7AA"> 

—14.87903 


500 


A OA _l_ A AO 

0.30 zb 0.02 


10.94 


19.6 


56.74 


1.09 


164.55 


25 


296.2156/ 


1/1 /ICC 

—14.84455 


500 


A A 1 I A AO 

0.41 zb 0.03 


12.76 


31.2 


77.53 


0.80 


241.65 


26 


on<i i t 1 n 
296. 13 132 


1 a ^noc c 
—14.69955 


500 


2.19 zb 0.1 / 


1 oo 
/.52 


oo o 
11 .1 


410.40 


0.56 


206.02 


27 


oa/; o/Coi£/i 

296.26364 


i a Anono 

—14.90800 


500 


O T£ 1 A 1 O 

2.36 zb 0.18 


6.94 


22.2 


449.03 


0.61 


173.41 


28 


296.2/594 


1 A Al A AA 

—14.92400 


500 


A OO I A AO 

0.23 ± 0.02 


10.50 


15.9 


43.85 


1.29 


138.64 


29 


oa^ ooco/; 

296.22586 


1 A T7£ n 

—14. 7763 7 


500 


A TO I A t\£L 

0.72 zb 0.06 


9.88 


24.8 


136.12 


0.78 


195.75 


30 


OA*C 0"7 1 1 O 

296.2/1 18 


1 yl OA 1 1 A 

— 14.89139 


500 


A AO I A A"7 

0.93 ± 0.07 


8.03 


18.6 


176.13 


0.84 


153.02 


3 1 


296.26489 


1/1 OO 1 /CO 

—14.8/163 


500 


OOO I A 1 1 

2.88 ± 0.21 


4.77 


1 1.6 


548.05 


0.81 


99.22 


32 


OA£ lO/iAc 

296.22665 


1 a ooooo 
— 14.5 / /2 / 


250 


1.15 ± 0. 1 / 


(~ O/l 

0.84 


O 1 A 

21.0 


A O A A A 


0.64 


165.29 


1 i 

3 J 


296.22006 


1/1 O/l £01 

— 14.540/z 


500 


3.11 ± O.Zj 


(~ /1A 

0.40 


22.0 


610.64 


0.57 


170.91 


34 


296.2271 1 


i /i Tn i c 

—14. /2315 


500 


1.56±0.12 


7.72 


22.3 


295.83 


0.68 


175.53 


35 


or* /; oaooa 

296.20230 


1/1 1"£1 Ol 

— 14./618/ 


250 


A in 1 A AO 

0.20 ±0.02 


12.45 


20.6 


37.25 


1.18 


173.04 


36 


296.19022 


-14.85646 


500 


1.49±0.12 


6.97 


17.8 


283.35 


0.77 


145.36 


37 


296.31396 


-14.97089 


500 


1.55±0.12 


6.52 


15.9 


294.40 


0.80 


131.66 


38 


296.20490 


-14.75966 


250 


1.37±0.10 


5.74 


11.6 


260.40 


0.97 


101.32 


39 


296.23492 


-14.72453 


95 


0.80 ±0.06 


7.49 


15.1 


152.41 


0.97 


128.27 


40 


296.23090 


-14.75460 


500 


1.10±0.09 


7.59 


18.2 


209.69 


0.82 


148.76 


41 


296.23541 


-14.73430 


500 


1.16±0.09 


6.20 


12.4 


219.38 


0.98 


107.88 


42 


296.19272 


-14.86216 


500 


1.44±0.11 


6.14 


13.6 


273.66 


0.89 


115.82 


43 


296.24420 


-14.73231 


500 


3.20 ±0.25 


4.89 


12.8 


607.40 


0.75 


108.19 


44 


296.31229 


-14.75217 


500 


1.31 ±0.10 


6.20 


13.2 


248.96 


0.92 


113.13 


45 


296.25198 


-14.88950 


500 


1.93 ±0.17 


6.60 


18.2 


366.74 


0.71 


146.81 


46 


296.24365 


-14.87206 


500 


1.52±0.13 


6.41 


15.2 


289.23 


0.83 


126.88 


47 


296.21973 


-14.88315 


250 


0.25 ±0.02 


9.78 


14.4 


47.81 


1.33 


127.30 



NOTE. — The assignment of an indicative maximum age (Col. 4) for each cluster is discussed in § 14.21 The 
estimation of the structural parameters of the clusters is described in detail in §|7] 



described as fractal) in scales starting from the largest giant 
molecular cloud down to individual clumps and cores. Within 
this scheme, the cluster population of NGC 6822, as pre- 
sented above, should represent one specific length-scale of a 
whole spectrum of stellar concentrations in this galaxy, and 
the detected clusters themselves are most probably members 
of larger structures. In this section we uncover the complete 
structural spectrum of stellar clustering and quantify the hi- 
erarchical distribution of the blue stars in NGC 6822 with 
r SS 500 Myr, as they are selected in §[3] 

8.1. Detection of stellar structures in NGC 6822 

In order to search for stellar concentrations of various 
length-scale s we apply again the NN density method, as de- 
scribed in § 14.11 for the 10th nearest neighbor (NN). The se- 
lection of the 10th NN is made for reasons of consistency with 
the detection of star clusters in the previous sections. Our tests 
showed that the application of the NN method for larger num- 



ber of neighbors (e.g., 30 or 50) would 'smooth' the derived 
structures and unify many of them into single larger ones, 
loosing the ability to detect any fine-structure in the spatial 
distribution of stars. For the identification of stellar structures 
of different densities we select different density thresholds in 
the application of the NN method. Specifically, apart from 
the 3cr detection applied in § [5] for the identification of the 
cluster population of the galaxy, we apply the method for the 
lower density thresholds of 1 a and 2a above the average back- 
ground density level to identify structures that correspond to 
lower stellar density enhancements. For the detection of con- 
centrations that correspond to high stellar density we apply 
the NN method with higher (4er and 5<r) density threshold 
above the background level. 

The stellar structures identified with the 10th NN density 
method are shown in the density contour map of Fig. [10] 
Isopleths of different colors indicate the structures derived 
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FIG. 10. — The density contour map constructed with the application of the 10th NN density method on our photometric data of NGC 6822. Isopleths within 
different density levels in cr, drawn with different colors, signify the corresponding identified stellar structures in the galaxy. All statistically important stellar 
structures are detected only in the main part of NGC 6822 disk. Therefore, the drawn contour map covers only this part of the galaxy instead of the whole 
observed field-of-view for reasons of clarity. 



with the application of the method for different detection den- 
sity thresholds. The cluster population of the 3cr detection 
is shown with green contours. This map is constrained to the 
main part of the galaxy, where all significant stellar concentra- 
tions, including i ts clusters, appear and which is known to host 
star formation dCannon et alJl2006T) . From the map of Fig.flOl 
it can be seen that the low density threshold detections reveal 
large stellar structures, while the application of the method 
for higher density thresholds give rise to smaller and more 
compact stellar concentrations, which are actually located in 
the larger structures of the galaxy. This combination of high- 
density enhancements that correspond to small stellar systems 
and aggregates with low-density concentrations that represent 
large structures and s tellar complexes is a clear signspot of 
hierarchical structure ( Elmegr eenl20 1 01) . 

8.2. Dendrograms 

An intuitive way to illustrate hierarchical structures is 
through the so-called dendrograms, introduced as "struc- 
tur e trees" for the a n alysis of molecul ar cloud structure 
by iHoulahan & Scal o (1992), refined by Rosol owskv et al.l 
(120081) . A dendrogram is constructed by cutting the im- 



age at different thresholds and identifying connected areas, 
while keeping track of the connection to "parent structures" 
(on a lower level) and "child structures" (on the next higher 
level, lying within the boundaries of its parent). A geometri- 
cally perfect hierarchy would be represented by a dendrogram 
where each parent branches out into the same number of chil- 
dren at each level. 

We construct the dendrogram of the stellar structures de- 
tected in NGC 6822 with the NN density method for the var- 
ious density thresholds considered, i.e., 1 to 5a above the 
background density. In this dendrogram, shown in Fig. [TU 
the structures found at each density level are represented not 
only by the 'leaves' that end at the particular level, but by all 
branches present at that level. For example, at the 3er level, 
there are 47 branches of the dendrogram (regardless whether 
they end here, continue to a higher level, or split into two 
or more branches), corresponding to the 47 detected stellar 
clusters. This dendrogram demonstrates that, while there are 
few centrally concentrated structures with a single peak, most 
structures split up into several substructures over at least three 
levels. The combination of this dendrogram with the contour 
map of Fig. [10] illustrates graphically the hierarchical spatial 
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FIG . 11 . — Dendrogram of the stell ar structures identified at different levels 
in the 10th NN density map (Fig. llOh illustrating the hierarchical behavior of 
the clustering of the blue MS stars. 

distribution of stars younger than 500 Myr in NGC 6822. 

In the lower detection density threshold of la (black iso- 
pleths in Fig. |T0T > there are few structures revealed, three of 
them being major stellar concentrations that qualify as large 
stellar complexes. These structures comprise a number of 
smaller multiple concentrations seen in the 2<r red isopleths 
of Fig. [10] which most probably correspond to the so-called 
stellar aggregates of the galaxy. The detected clusters (green 
3cr isopleths) are actually members of these aggregates, ful- 
filling the typical image of hierarchical structuring of stars in 
a galaxy-scale. Higher density (4 and 5a) detections corre- 
spond to the condensed stellar density peaks, which are seen 
to be the most compact centers of the larger structures. 

8.3. Quantification of hierarchy 

Another popular method in studies of hierarchical cluster- 
ing is the Minimum Spanning Tree (MST) method, which 
may be applied with different cutoff lengt hs for the detec - 
tion of different length-scales (e.g., Bas tian et al.l20 07. 2009). 
The MST method is excellent in quantifying the hierarchy in 
the disco vered structures through the so-called Q parameter 
(see e .g., ICartwright & Whitworthll2004t ISchmeia & Klessenl 
2006). It allows to distinguish between clusters with a central 
density concentration and hierarchical stellar concentrations 
with possible fractal substructure. Large Q values (Q > 0.8) 
describe centrally condensed clusters having a volume density 
n(r) oc r~ a , while small Q values (Q < 0.8) indicate concen- 
trations with fractal substructure. Q is correlated with the ra- 
dial density exponent a for Q > 0.8 and anti-correlated with 
the fractal dimension D for Q < 0.8. We determine Q for all 
structures detected at all levels with more than 45 members, 
since for poorer concentrations the method becomes unreli- 
able due to increasing errors. Recent simulations of artifi- 
cial clusters have shown that for clusters with N <45 mem- 
bers, the expected error in the determin ation of the Q pa- 
rameter is of the order of less than 0.05 (Schmeja 201(J). Q 
varies s ystematically for structures showing a more elongate d 
shape (ICartwright & Whitworthl 120091; iBastian et al] 12 009), 
therefo re we apply the correction suggested by Bastia net al.l 
(2009) when necessary. 

The Q values for the structures detected at 1, 2, and 3a 
density thresholds are all but one in the hierarchical regime 
(Q < 0.8), reaching Q values as low as Q = 0.55, corre- 
sponding to a fractal dimension of D » 1.8. The detections 
at higher density thresholds did not provide stellar structures 
rich enough to allow such an analysis. The Q values for the 
five lcr structures with more than 45 stellar members range 



from 0.56 to 0.73, and those for the ten 2a structures from 
0.55 to 0.83. There are only two clusters detected in the 3a 
density threshold that include more than 45 stars, clusters No. 
1 and 2 in Table [JJ They are also found to be hierarchically 
structured with values of Q = 0.62 and 0.77 respectively. The 
average Q value for the structures revealed in all three density 
levels give evidence of fractal concentrations with Q = 0.70. 
This value demonstrates the hierarchical nature of the cluster- 
ing in NGC 6822. All but one of the analyzed structures are 
classified as hierarchical, i.e., they have significant (possible 
fractal) substructure rather than a smooth density gradient. It 
should be noted that the numbers of stars in the objects re- 
vealed at density levels with thresholds > 3(7 are too low to 
permit a meaningful Q analysis, but taking as example the 
two star clusters, for which this analysis was possible, most - 
if not all - of the 3 a clusters should be expected in the hierar- 
chical regime as well. 

Nevertheless, the unusual distribution of blue stars in the 
central part of NGC 6822 and their apparent hierarchy may 
be the result of differential extinction. For example, the two 
vacant areas to the east and west of the central S-shaped dis- 
tribution of blue stars (Fig.fTOl. may be caused by high levels 
of extinction. Indeed, a comparison of the spatial distribution 
of the cool component of neutral hydrogen with that of the 
blue stars shows an anticorrelation on large scales (dBW06, 
their Fig. 7). However, while this extinction seems to af- 
fect the large-scale stellar distribution (at la density level), 
giving it the S-shape, it does not seem to affect the hierarchi- 
cal distribution of smaller stellar structures detected in higher 
density levels (2 and 3a). Concerning the smallest detected 
structures (at 4 and 5a), the resolution of the available obser- 
vations of the ISM in NGC 6822 does not allow the detection 
of any patchy extinction in smaller scales. As such, we can- 
not assess how small-scale differential extinction may affect 
the apparent hierarchy of stellar structures, and therefore we 
certainly cannot exclude it as a possible bias. 

In this case the derived Q values would indicate that the 
stars were in fact not in total hierarchically spread. It is also 
worth noting that if the largest stellar structures, shown in 
Fig. [TTJ are approximately 2D (i.e., much longer and wider 
than deep) the Q values returned by the MST method should 
be interpreted differently, since the crossover between 'frac- 
tal' a nd 'centrally concen trated' in the 2D case becomes 0.72 
(see, [gastian et al. 2009), and thus the structures may be 
found to be less hierarchical. 

8.4. Correlation of number of stars and radii 

The correlation between the number of stars and radius seen 
for the clusters of NGC 6822 detected at the 3a level (§ r7~TT i 
is also observed for the structures we identify with the ap- 
plication of the NN method with different detection levels. 
This correlation is shown in Fig.[T_2]for detections with den- 
sity thresholds of 1, 2, 4 and 5er. The number of stars cor- 
relates with the radii of the structures with about the same 
slope of k « 0.53. Only the 5cr structures show a differ- 
ent beha vior, but this plot su ffers from the small number of 
objects. Adam s et al.l (|2006) found a slope of k = 0.54 for 
young Galactic clusters, very close to the slope we find for 
the detected ste llar structur e s and clusters in NGC 6822. As 
pointed out by Alle n et all (T2007), these results are similar 
to that of Lar son's relations for molecular clouds, p ~ R~ l 
dLarsonl 19811) . implying a constant column density of the gas 
in molecular clouds. 
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FIG. 12. — Correlation between the number of detected stars and r equ i v for the stellar structures of NGC 6822 identified with the NN method for various density 
detection lev els. The corresponding slopes, k, and correlation coefficient s, q are also given. This correlation is similar to that we find for cluster population of 
the galaxy (§ 17. FJ and consistent with findings for Galactic clusters (see § 17.11 . 



8.5. Current Star formation in NGC 6822 

ICannon et alj (120061) performed Spitzer Space Telescope 
imaging of NGC 6822 and combined it with Ha, Hi and ra- 
dio continuum observations to study the nature of the emis- 
sion in these wavebands on spatial scales of ~ 130 pc. They 
found strong variations in the relative ratios of Ha and IR flux 
throughout the central disk of the galaxy, and that the local- 
ized ratios of dust to Hi gas are about a factor of 5 higher than 
the global value of NGC 6822. These authors identified 16 
IR emission complexes in the central part of the galaxy, six 
of which correspond to detections in the radio continuum im- 
age. The major Hll regions, i.e., the strongest Ha sources of 
the galaxy are also luminous in the IR and identified among 
the IR emission complexes, but in general Ha and IR emis- 
sion are not always cospatial. The Hi distribution in the same 
area contains mostly high surface brightness (column densi- 
ties <; 10 21 cm -2 ), and it is clumpy with dense comps of 
neutral gas surrounding various emission peaks in other wave- 
bands. 

It is wort hwhile to compare the star-forming complexes 
identified bv lCannon et alj (120061) with the structures revealed 
from our study. The IR emission complexes are indicated by 
circles in Fig. [13] overlaid on our NN density map of the cen- 
tral part of NGC 6822 disk. This map appears more smooth 
than those shown before because we constructed it by apply- 



ing the NN density method for the more gross number of the 
50th NN, in order to achieve an angular resolution compara- 
ble to that provided by the nebular and dust emission obser- 
vations. The comparison of our stellar structures with the IR- 
bright regions of Cannon et al. on the map of Fig. [T3l shows 
a correlation of the stellar density peaks with the maximum 
IR emission only in the northern part of the galaxy, in partic- 
ular, their regions 9 and 10 coincide with peaks in the stellar 
density. In the southern part they are significantly displaced 
from each other. Here, current star formation takes place at 
the western rim of the region of high stellar density. This is 
reminiscent of the spatial distribution of stellar birth in spiral 
arms or tidal tails, and points to wards an external interaction 
(see also Karampelas et al. 2009). On the other hand, this spa- 
tial disagreement between peaks in IR emission and in stellar 
density in the southern central region of NGC 6822 may sim- 
ply suggest that young stars are not yet revealed in their natal 
clusters, and we only detect the most evolved stellar concen- 
trations in this part of the galaxy. 

9. DISCUSSION AND CONCLUSIVE REMARKS 

The formation of stars and star clusters in the interstel- 
lar matter (ISM) is controlled by the complex dynami- 
cal interplay between self- gravity and supersonic turbulence 
dMac Low & Kl essen 2004). On large scales, the turbulent 
motions in the ISM are highly supersonic with Mach numbers 
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FIG. 13. — The 50th NN density map of the blue MS stars in the central re- 
gion of NGC 6822 (in logarithmic scale) with blue circles overlaid indicating 
the IR-bright complexes identified by Cannon et al. 12006) with Spitzer. 

up to several tens in giant molecular clouds. On scales of in- 
dividual star-forming cloud cores, the turbulent velocity field 
becomes sub sonic, with a w ell-defined power-law connecting 
these scales (lLarson|["l981l) which could be the result of the 
scale-free nature of the turbulen t cascade (Elmegr een & Scalol 
120041: IScalo & E lmegreen 2004). The resulting morphologi- 
cal structure has often been interpreted in terms of hierarchi- 
cal fr actals (e. g.. lElmegreen & Falgarone||1996t IStutzki et alJ 
1998; see also lFederrath et al II2009L l20ld> . 

Turbulence plays a dual role. On global scales it pro- 
vides support, while at the same time it can promote col- 
lapse locally. It creates strong density fluctuations with 
gravity taking over in the densest and most massive re- 
gions, where co llapse sets in to build up individual stars 
dKlessenll2001al) . Together with the thermodynamic prop- 
erties of the gas and magnetic fields this regulates the 
fragmentation behavior in star forming clouds. The prop- 
erties of young star clusters are thus directly related to 
the statistical characteristi cs of the underlying turbulenc e 
in the star-forming gas (iBallesteros-P aredes et al. 2007). 
This holds for the stellar IMF dKlessen & Burkerj 120001 
I200U iHennebelle & Chabrierl 120081 12009b. the proto stellar 
accretion rates (IKlessenll2001bt ISchmeia & Klessenl [2 004). 
the star-formation efficienc y and timescale dKlessen et alJ 
20001 iHeitsch et all [200ll IVazquez-Semadeni et all 120031; 
Krumholz & McKee 2005), or the question of when and 
where dense clusters with high-mass stars form and when and 
where to e xpect a more distributed popula tion of lower-mass 
stars (e. g.. IVaz quez-Semadeni et al. 20091). 

Quite a number of different measures have been pro- 
posed to support and statistically characterize the link be- 



tween the morphological and kinematical structure of the 
ISM and the population of star clusters that form in it. 
lEfremov"& Elmegreen ( 1998) find a correlation between the 
separation of stars and star clusters and their age, and inter- 
pret that result in term s of the linewi dth-size relation in the 
ISM first discussed by Larsonj d 198 lb . As velocity and size 
converts into a timescale, they conclude that the formation of 
molecular clouds and stellar birth in their int erior proceeds on 
timescales of a few turbulent crossi ng times dElmegreen et alJ 
2000; Elmegreen 2007) (see also Ballesteros-Paredes et al. 
19991: lHartmann et alJ l2001h iGlover et all 120101) which~~is 
a signpost of star formation controlled by interstellar tur- 
bulen ce (lElmegreen & Efremovl 119971 ; Mac Lo w & Klessenl 
2004). Further evidence is provided by the fact that the mass 
spectrum of molecular clouds and clumps within the clouds 
follows a power law, dN/dM oc M~ a , with an exponent in the 
range 1.5 < a < 2.0 (Stutzki & Guestenll990t I Williams et alJ 
1994; KrameretaDll998j) which" is very similar to the mass 
spectra of young star clusters as discussed in § [6j The value 
of a = 1.5 ±0.7 that we derive for the young clusters in 
NGC 6822 is fully consistent with this picture. 

As the ISM exhibits a very complex, hierarchical morpho- 
logical structure, it stands to reason that stars follow a similar 
spatial pattern. We have studied that aspect in §[8]in terms of 
the distributions of blue stars in individual stellar structures. 
Using the 10th nearest neighbor map to identify clusters, we 
find that indeed these stars in NGC 6822 can be grouped into 
larger and larger aggregates in a hierarchical fashion when we 
vary the detection th reshold. Similar re s ults are found in the 
Magella nic Clouds (iGieles et al.1 120081: iBastian et alJj2009b 
see a lso iHunteretal J 120031) . M33 dBastian et al J 120071) . and 
M51 dBastian et al.ll20 05). We note that the Q-parameters de- 
rived for the individual star clusters in NGC 6822 indicate that 
the stars in the clusters themselves are hierarchically struc- 
tured. It is a common feature of young star clusters that they 
reveal a high degree of sub-clustering when observed with 
sufficient resolution and sensitivity t o identify individual stars 
down to the peak of the IMF (e.g.. |Car twright & Whitworthl 
2004, 2008l l2009tlSchmeia et alj2008ll2009t) with values that 
are con sistent with numerical calc ulations of star-cluster for- 
mation (Schmeia & Klessen 2004). We conclude that we find 
direct evidence, that the blue stars in NGC 6822 exhibit a hi- 
erarchical spatial pattern from the scales of individual objects 
all the way up to the scale of the galaxy as a whole. 

It is interesting to note in this context that when it 
is indeed true that the spatial distribution of stars and 
star clusters traces the statistical properties of ISM tur- 
bulence, then the absence of a clear break in the hi- 
erarchy could indicate that turbulence is driven on kpc- 
scales. Possib le mechanism s are global gravitational in- 
stability (e.g.. iLi et al.1 12005b . the magnet o-rotational insta- 
bility (MRI) in the disk of NGC 6822 dBeck et ail [19961; 
iBalbus & Hawlevlll998l;|Piontek & Ostrikerll2007l) . or the ac- 
cretion of fresh gas t hrough the halo in form of cold streams 
or gas-rich satellites dSantillan et alj2007fclSancisi et alj2008l; 
iKlessen & H ennebella 120101) . Indeed, the "Northwestern 
Cloud" clearly visible in Fig. Q] has been speculated to 
be a separate system, possibl y interacting with NGC 6822 
dde Blok & Walterl2000l. 120031) . This interaction could also be 
responsible for a pos sible increase of the star f ormation rate 
in the past 100 Myr dHodgd ll980^ see however G allart et ail 
1996a, for a different opinion). 

10. SUMMARY 
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In this paper we present a thorough investigation of the 
cluster population with age ^ 500 Myr and the hierarchy 
in the spatial distribution of main-sequence stars of this age 
in the Local Group dwarf irregular galaxy NGC 6822. Our 
observational material comprises optical imaging from the 
8.2-m SUBARU Telescope, providing the most complete point 
source catalog of the galaxy in terms of dynamic range and 
spatial coverage. 

Star clusters are identified with the application of the 
Nearest-Neighbor (NN) density method for the 10th NN on 
the blue main-sequence stars in NGC 6822. Fourty-seven dis- 
tinct concentrations with > 3 stellar members that are found 
with our method to have 10th NN density values 3a above 
the average background density are identified as the star clus- 
ters of the galaxy. The physical dimensions, stellar density 
and limiting radii of the clusters are defined by this density 
limit. The size distribution of the detected clusters can be very 
well approximated by a Gaussian peaked at ~ 68 pc (with 
<r~18pc). 

The total stellar masses of the clusters are estimated by ex- 
trapolation of the total observed stellar mass function of all 
clusters to the lower stellar masses, assuming that it has the 
shape of the average Galactic mass function ( Kroupa| [2002l) 
down to 0.1 Mq. The derived cluster masses resign in the 
range 10 3 - 10 4 M©. Their distribution follows very well 
a power-law with index ~ -1.5 ±0.7, somewhat shallower 
but consistent with the Cluster Initial Mass Function and the 
mass spectra found for clusters in other Local Group galax- 
ies. Structural parameters, such as their Spitzer radius, the 
correspondent 3-dimensional stellar density and velocity dis- 
persion, as well as their crossing and the two-body relaxation 
times are computed for the clusters from their total masses 
and stellar numbers. We find a correlation between the radii 
and stellar numbers of the detected clusters, which follows a 
power law similar to that found for Galactic embedded clus- 
ters. 

We study the clustering behavior of the blue stars in the 
considered age-limit with the application of the NN density 
method and the detection of stellar structures in different den- 
sity levels, other than 3a used for the establishment of the 
clusters. We identify, thus, stellar concentrations of various 
length-scales, clusters being the representative of one of them. 
High-density stellar structures are revealed from the applica- 
tion of the NN method with higher (4a and 5a) density thresh- 
old above the average density level, while the larger loose stel- 
lar structures, which include the denser ones, are recognized 
by the application of the method for lower (la and 2a) den- 
sity thresholds. All significant stellar concentrations, includ- 
ing the clusters, appear in the main (star-forming) part of the 
galaxy. The density maps constructed with this technique pro- 
vide clear indications of hierarchically structured stellar con- 
centrations in NGC 6822, in the sense that small dense stellar 
concentrations are located inside larger and looser ones. We 
illustrate this hierarchy by the so-called dendrogram, or struc- 
ture tree of the detected stellar structures in NGC 6822, which 
demonstrates that most of the detected structures split up into 
several substructures over at least three levels. 

The majority of the detected star clusters in NGC 6822 
has ages between 250 and 500 Myr, which is longer than the 
crossing time of the galaxy, the timescale on which hierarchi- 
cal cluster dis tributions in the LMC and SMC are erased (e.g., 
Bastian et al. 2009). As a consequence, the finding that these 



clusters have preserved their structure is in itself very inter- 
esting. However, our catalog of "stellar clusters" comprises 
large stellar groups with sizes larger than the typical pc-scale 
Galactic clusters; the smallest cluster in our catalog is about 
20 pc in size. As such, their large sizes is a possible explana- 
tion why these clusters have preserved their structure longer 
than one crossing time of the galaxy. 

We quantify the hierarchy of these structures with the use 
of the Minimum Spanning Tree method and the Q param- 
eter. We find that the Q values for the structures detected 
at 1, 2, and 3a density thresholds reside in the hierarchical 
regime (Q < 0.8), reaching Q values that correspond to a 
fractal dimension of D ps 1.8. The correlation between the 
number of stars and radius seen for the clusters of NGC 6822 
is also observed for the larger stellar structures, identified at 
the la and 2a density levels, as well as for the smaller dense 
concentration identified at the 4a density level. It is worth 
noting that some of the large high-density stellar concentra- 
tions, particularly in the northern part of the star-forming por- 
tion of the galaxy, coincide with IR-bright structures earlier 
identified with Spitzer, associated with high column density 
( 10 21 cm" 2 ) neutral gas, indicating structures that currently 
form stars. 

The clusters of NGC 6822, which we identify at the 3a den- 
sity level of the NN method, are found to be themselves inter- 
nally structured. Considering the ages of these clusters, this 
result is also interesting, because as found for small compact 
clusters, any substructure is expe cted to have been rapidly 
erased on very short time-scales ( I Allison et a l. 2010). How- 
ever, the clusters we detect in NGC 6822 are large stellar 
groups, which may have survived strong disruptions for peri- 
ods of time longer than their c rossing times. Previous simula- 
tions of fractal star clusters by Goodwin & Whitworth (2004) 
have shown that even an initially homogeneous cluster can 
develop substructure, if it is born with coherent velocity dis- 
persion. As a consequence, the substructure we observe in 
our clusters may not be the present exhibition of their original 
'primordial' substructure, but one induced later. The same 
authors find that the velocity dispersion is a key parameter, 
determining the survival of substructure, which can last for 
several crossing times. Therefore, it is possible that at least 
some of our clusters may have sustained their initial substruc- 
ture. In favor to this, spatial substructure has been identi- 
fied in open clusters of t he Milky Way as old as ~ 100 Myr 
dSanchez & Alfaroll2009l) . 

The morphological structure of the stellar concentrations 
we identified in NGC 6822 resembles the structural behavior 
of the interstellar matter, which in principal is hierarchical and 
dominated by turbulent motions. As such, this coincidence 
between the clustering of stars and the hierarchical structure 
of the ISM, as it is demonstrated in the present study, suggests 
that turbulence may play the major role in regulating clustered 
star formation in NGC 6822 on pc- and kpc-scales. 
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